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We present cluster-DMFT (CTQMC) calculations based on a downfolded tight-binding model 
in order to study the electronic structure of vanadium dioxide (VO2) both in the low-temperature 

Q (-^1) cind high-temperature (rutile) phases. Motivated by the recent efforts directed towards tun- 

ing the physical properties of VO2 by depositing films on different supporting surfaces of different 
orientations we performed calculations for different geometries for both phases. In order to investi- 

G gate the effects of the different growing geometries we applied both contraction and expansion for 

jH the lattice parameter along the rutile c-axis in the 3-dimensional translationally invariant systems 

O 

Q miming the real situation. Our main focus is to identify the mechanisms governing the forma- 

tion of the gap characterizing the M\ phase and its dependence on strain. We found that the 

vQ increase of the band-width with compression along the axis corresponding to the rutile c-axis is 

00 

QQ more important than the Peierls bonding-antibonding splitting. 
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PACS numbers: 71.15.Ap, 71.27.+a, 71.30. +h, 79.60.-i 



I. INTRODUCTION 

Vanadium oxides compounds exhibiting exotic transport phenomena are subjects of ex- 
tensive interest. In particular vanadium dioxide, VO2, undergoes a first-order transition 
from a high-temperature metaUic phase to a low-temperature insulating phase at almost the 
room temperature (T = 340 K) [T]. There are intensive efforts around the world to make 
devices such as switches, transistors, detectors, varistors, phase change memory, exploiting 
the unique properties of VO2. 

At low-temperature VO2 has a simple monoclinic (Mi) structure with space group P2i/c 
(Ml phase) while at high temperature it has a simple tetragonal lattice with space group 
P4:2/mnm rutile (i?-phase) as displayed in Fig. [1} The lattice structures of the two phases 
are closely related as emphasized in Fig. [T]by showing similar wedges of the different lattices: 
the Ml unit cell is similar to the R unit cell, when the latter is doubled along the rutile c-axis. 
For the sake of simplicity, we will use the notation c-axis for both the Mi and R phase when 
referring to the axis equivalent to the rutile c-axis (in Mi phase this axis is sometimes called 
a-axis). The Mi phase is characterized by a dimerization of the vanadium atoms into pairs, 
as well as a tilting of these pairs with respect to the c-axis [21 [3] as indicated by showing 
bonds between the paired V atoms in Fig. [1]d. 

The electronic and transport properties are dramatically different for the two phases. 
The resistivity jumps by several orders of magnitude through the phase transition, and the 
crystal structure changes from i?-phase at high-temperature to monoclinic Mi-phase at low- 
temperature [H H] . While the rutile phase is a conductor, the Mi phase is an insulator with 
a gap of ~ 0.6 eV [5] at the Fermi energy. VO2 has also attracted a great deal of attention 
for its ultrafast optical response, switching between the R and the Mi phases [BHS]. Despite 
the large number of experimental [5l [814T8] and theoretical [T9H27J studies focusing on this 
material the physics driving this phase transition and the resulting optical properties are 
still not identified undoubtedly. In the theoretical works using LDA [U |211 ES] the formation 
of the gap was not found. The single site DMFT approach [261 EZj is known to correctly 
describe the metallic phase but can not take into account the formation of the bonding 
states of the dimers in the insulating phase which requires the cluster DMFT method [19] . 
As shown in Ref. [21] the fully interacting one-electron spectrum of the Mi phase can be 
reproduced by an effective band structure description. The optical conductivities calculated 




FIG. 1: (Color online) a): The high temperature metallic rutile structure, b): The low temperature 
insulating Mi structure. Large (red) spheres: V atoms, small (blue) spheres: O atoms. Thin black 
box in the rutile structure emphasizes the similarity between the two phases. Thick black boxes: 
the unstressed unit cells, thick gray (orange) boxes: stressed unit cell. Arrows: the directions of 
lattice vectors (note the names of the axes differ from the usual notation [24J). The heights of the 
unstressed (stressed) unit cells are cq (c = rco) for the rutile and cq (c = rco) for the Mi structure. 
Green lines connect the dimerized V atoms. The actual stress in theory and experiments is 10 
times smaller than depicted in the figure. 



by using cluster DMFT shows good agreement with the experiment for both phases [43]. 

The metal-insulator transition (MIT) of VO2 is usually attributed to two different physical 
mechanisms. One of them is the Peierls physics, i.e. the dimerization of the V atoms along 
the rutile c-axis [28], and consequently opening of the gap in the Brillouin-zone reduced 
in the direction along the c-axis. In this case, the gap opening can be explained in the 
framework of effective one-electron theories. The other one is the Mott mechanism where 
the gap opens due to the strong Coulomb repulsion between the localized V-Sd orbitals 
[29] and the related dynamical effects. Understanding in detail the interplay and relative 
importance of both the Peierls and the Mott mechanisms for the electronic structure is 
crucial for controlling this material with an eye towards applications. For example, whether 
the driving force of this transition is electronic (i.e. occurring on femtosecond timescales) 



or structural (occurring on the picosecond timescale) is important to understand the speed 
of the switching from the Mi to the rutile phase. From the perspective of apphcations, in 
order to control the properties of vanadium dioxide, it is essential to identify the effects of 
compressive and tensile stress resulting from the various substrates, on which the films are 
deposited p[TO]. 

The experimental results of Ref. [9J showed that compressive uniaxial strain (along the 
c-axis) stabilizes the metallic phase. This result cannot be explained by applying a simple 
Peierls picture exclusively. The Peierls mechanism predicts that a compression along the 
c-axis increases the splitting between the bonding and antibonding state, formed by the 
combinations of particular 3d orbitals residing on the different V atoms in a dimer, and 
hence would promote insulating state. The Peierls mechanism thus increases the tendency 
to open the gap in the Mi phase, which would stabilize the insulating phase at the expense of 
the metallic phase. In this picture, the transition temperature would increase under uniaxial 
compressive strain, opposite to experiment [H]. 

The LDA calculations fail to reproduce the gap in the Mi phase [21]. This suggests 
that the correlations are important in this material. Motivated by this and the above 
mentioned experiment we have examined the influence of strain on the electronic structure 
of VO2 microscopically applying the dynamical mean-field theory j30] combined with the 
local density approximation of density functional theory (LDA+DMFT) [31]. We carried 
out LDA+DMFT calculations of VO2 under strain. We established that in addition to the 
Peierls distortion, other factors like the position of the e^ band and the band-width of the 
bonding state in the Mi phase showing definite sensitivity to the strain, play a defining 
role in driving the MIT in VO2. Theoretical predictions for the strain dependence of many 
spectroscopic quantities are made, including the photoemission, the optical conductivity, 
the inverse photoemission and XAS in the R and Mi phases. The insights achieved in this 
study together with the computational machinery developed, will serve as a basis for rational 
material design of VO2 based applications. 

II. METHOD 

The unit cell of the rutile structure contains two vanadium and four oxygen atoms while 
the Ml structure contains four vanadium (2 dimers) and eight oxygen atoms. The calcu- 



lations in the rutile phase were performed with the doubled unit cells (four V and eight 
O atoms) to allow the formation of a bonding state between the V atoms separated along 
the c-axis, a mechanism which has dramatic effect in the Mi phase. The lattice structure 
parameters for the rutile (referred as oq, &o and cq in the following) and Mi cases (oq, &o, 
Co and /3o) were published by McWhan et al. in Ref. [37] and Kierkegaard and Longo in 
Ref. p], respectively. Note that our notation of the axis differs from the published data [2^ . 
in order to emphasize the similarity of the R and Mi structures. 

The strain was applied by changing the lattice constants but not the internal parameters 
of the atomic positions. The strain is characterized by the ratio between the experimental 
lattice constant along the c-axis of the unstrained crystal, and the one used in the calculation, 
r = c/cq in the rutile and r = c/cq in the Mi case. The applied strains correspond to 
r = 0.98, 1.00 and 1.02 ratios both for the R and Mi phases in order to make calculations 
which are compatible with the experiment in Ref. [9]. The lattice parameters were changed 
with the constraint to preserve the volume of the original unit cell and to form a structure 
to closely resemble the original structure. This is illustrated in Fig. [T]for both phases. The 
constant volume constraint in the rutile phase requires for the other two lattice constants 
a = ao/\/r and b = 6o/^A'• For the Mi structure, two other constraints were introduced: 
the two ratios {bsm/3)/a = 0.9988 and — 2(6sin/3)/c = 1.0096 were kept constant. Hence we 
used the following lattice constants and angles c = rco, a = So/a/t, tan(/3) = tan(/3o)/r^/^ 
and b = boy^r"^ cos'^{j3q) + sin^(/3o)/r. 

For the i?-phase, we used a local coordinate system introduced by Eyert in Ref [211 as 
shown in Fig. [2J In the rutile geometry the x-axis is parallel to the rutile c-axis while 
the z-axis points along the [110] direction, pointing to the adjacent oxygen atom in the 
vanadium plane. In the Mi phase, our local coordinate system is somewhat tilted to align to 
monoclinic geometry. The symmetry classification of the electronic orbitals is also adopted, 
namely d^z^_r2 and d^y stand for the Cg states and d^z, dyz and d^2_y2 for the t2g orbitals. 
The t2g manifold is divided further: the notation aig is used for d^^-yi state and e^ for the 
dxz and dyz states. In this choice of the local coordinate system the dimerization strongly 
affects the d^^-yi states, which have large electron density along the V-V bond. 

First, a self-consistent LDA calculation was performed using the linear muffin-tin orbitals 
method combined with the atomic sphere approximation (LMTO-ASA) [32] • For a satis- 
factory description of the interstitial region 48 empty spheres were added in the LDA-ASA 
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FIG. 2: (Color online) The sketch of vanadium t2g orbitals {dx2_y2, dxz, dyz) together with the 
applied coordinate system and the V-V pairs considered in the rutile phase calculations. Large 
(red) spheres: V atoms, small (blue) spheres: O atoms. 

calculations. In the next step we determined a downfolded parameters of the low energy, 
effective Hamiltonian, H^jf, including only the 3d t2g subset of electronic states lying in the 
vicinity of the Fermi energy. The effective Hamiltonian is constructed as: 






(1) 



with off-diagonal hopping (t) and diagonal one-electron energy (e) parameters belonging to 
different sites {i) and states (a = d^z, dyz, d^^-y^) with a spin character a. The applicability 
and accuracy of the above downfolding method is determined by the mutual positions and 
characters of the electronic bands. In the studied case of VO2 both in the Mi and R phases 
the bands in the proximity of the Fermi level have mainly t2g character and are well separated 
from both the low-lying 2p-bands of the oxygen, and from the Cg bands of vanadium, due 
to the strong crystal-field splitting (~ 3.0 — 3.5 eV [2l])- The crystal field splitting is due 
to oxygen octahedron surrounding the individual V atoms. This mapping of the LDA band 
structure onto a tight-binding type Hamiltonian provides a reasonably good description 
of the electronic structure as compared to the LDA serving as a good starting point for 
the following DMFT calculations. The evolution of the intra dimer V-V hopping matrix 
elements corresponding to the three different t2g orbitals induced by the strain are shown 
in Fig. |3j Although, the relative changes in the parameters are small, the manifestations of 
these changes in the low energy spectral function, and therefore in the transport properties, 




FIG. 3: (Color online) The intra dimer V-V hoppings corresponding to the t2g orbitals shown in 
Fig.^as a function of strain. The left (right) panel shows the parameters in the rutile (Mi) phase. 
The direction of strain is depicted in Fig. [T] 

are significant. 

The effects of the electron-electron interaction beyond LDA are treated by means of the 
dynamical mean-field theory [SUl EI] • The coulomb interaction was taken into account for the 
V-V dimers by using the cluster extension of the DMFT [3T]. In this approach the quantities 
describing the electronic states of a cluster (Green's function, self-energy, etc.) are matrices 
of both the site and electronic state indices. Having two sites (two V atoms) and three states 
{t2g states) the resulted matrices have dimensions of 6 x 6 for non-magnetic calculations. For 
solving the impurity problem, we used the continuous time quantum Monte Carlo method 
(CTQMC) [23] at the temperatures T = 232 K and T = 390 K, below and above the critical 
temperature, for the Mi and R cases, respectively. In the present calculations we assumed 
on-site Coulomb interaction on the V atoms written as: 

i=l,2 aa' aa' 

excluding the terms with similar orbital and spin character at same time. We fixed the 
parameter U at U = 2.2 eV, to reproduce the experimentally measured gap in the Mi 
phase, which shows a weak temperature dependence between 100 and 340 K, varying between 
~ 0.75 and ~ 0.6 eV [5l |3l]. The strong dependence of the gap size on the U parameter is 
shown in Fig. |4j 

The linear U dependence of the renormalization factor Z calculated from the real part 
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FIG. 4: (Color online) The variation of the t2g density of states of the Mi phase as a function of 
the applied Coulomb repulsion parameter, U. Note the high sensitivity of the gap on U. 



of the self-energy of the different orbitals (a) as Za 



dReT,a{(jj)/dijj) ^ and the elec- 



The Z value 
19] . An alter- 



tronic specific heat, 7 = ^ pa{0)/Za, of the rutile phase are shown in Fig 

a 

obtained at U = 2.2 eV is in good agreement with the one published in Ref? 
native method to determine the value of Z factor is to take the ratio of the experimentally 
measured plasma frequency in the rutile phase {Up^^ = 2.75 eV)[Tl] and devide it by the 
band theoretical LDA-LAPW calculation [35] (w^^^ ^ 4.1 eV), i.e., u^^p/loI^^^ ^ 0.67, 
which agrees well with Z ^ 0.62 in the case of f/ = 2.2 eV. Although the U was chosen to 
reproduce the gap in the Mi-phase, it is satisfying that this value is also compatible with 
the alternative estimation, based exclusively on the optical and thermodynamic properties 
of the rutile phase. 

In order to achieve a structured (almost diagonal) self-energy matrix for the cluster, we 
used a basis of symmetric (s) and anti-symmetric (as) combination of the states localized 
on the individual V atoms of the dimers defined as 
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To obtain the physical properties at real energies we performed analytic continuation to the 
real axis using a recently developed method of expansion in terms of modified Gaussians 
and a polynomial fit at low frequencies, as described in details in Ref. [36] . 
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FIG. 5: (Color online) Renormalization constants Zq, = (1 — dKeTja{uj) / dco) ^ of the different 
orbitals in the R phase as a function of U. Inset: electronic specific heat, 7 = ^Pq(0)/Zq, 

a 

calculated at different temperature and values of U. 
III. RESULTS AND DISCUSSION 



Fig. p^ displays the LDA t2g densities of states (DOS) of the V atoms in both the rutile 
and Ml phases for different r ratios (r = 0.98, 1.00, 1.02). The total band-width (~ 2.6 eV 
for both phases) and also the fine details of the DOS obtained by the LDA-LMTO agree 
well with previous studies [211 EH] • The minor discrepancies are probably a consequence of 
the different electronic structure method or the slightly different geometry. 

The effect of the stress applied along the rutile c-axis on the electronic DOS can be 
clearly seen. The position of the aig peak is changed considerably, while the total band- 
width is only slightly changed (< 0.1 eV) under application of stress. Only the d^i-y-i 
orbital changes significantly, while the two e^ orbitals remain mostly unchanged. This is 
due to the sensitivity of the overlap integral of (i^.2_y2 states along the c-axis. It is clear 
from Fig. M that the splitting of the d^^-y^ peaks strongly increases with decreasing lattice 
parameter c roughly following the linear dependence of the intra-dimer hopping parameter 
t^2_y2^^2_y2 on ratio r, shown in Fig. [3] However, even under compressive stress of r = 0.98, 
the splitting is not large enough to open a gap in the Mi phase. The value of the intra- 
dimer hopping parameters corresponding to the dx2-y2 states, which plays a significant role 
in the formation of the electronic structure of VO2, are —0.30 eV and —0.61 eV for the R 



and Ml cases, respectively. These values are in a good agreement with published values 
of Ref. [ini One can observe that the splitting of the d^^-yi peaks are larger for the Mi 
geometry (~ 1.4 eV) than for the rutile structure (~ 0.76 eV). This can be attributed to 
the reduced V-V distance along the c-axis due to the dimerization. The splitting of the 
bands can be roughly approximated by 2tj.2_j^2 ,j.2_y2 for the different r ratios. This behavior 
resembles the bonding and anti-bonding splitting of a dimer molecule, suggesting that the 
splitting of these states is determined mainly by the intra dimer hopping, especially in the 
Ml case. For the e^ states, this correspondence is less clear, showing the importance of the 
inter-dimer hoppings. It is worth to note that the LDA calculations can capture only the 
Peierls physics, which alone is not sufficient to explain the formation of the gap even if a 
reduced V-V distance is considered. 
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FIG. 6: (Color online) LDA partial DOS of the t2g states of a V atom in the rutile (left) and Mi 
(right) phases in the proximity of the Fermi energy with different r ratios: r = 0.98 (dotted, red), 
1.00 (full, black) and 1.02 (dashed, green). Note the large splitting of the dx2_y2 states and the 
barley changed DOS at the Fermi level. 

To trace the effect of the electron-electron interactions e.g., appearance of the gap in the 
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FIG. 7: (Color online) Comparison of Y-3d, 0-2p and the total DOS of LDA+DMFT calculation 
in rutile (left) and Mi (right) phases (r = 1.00) with angle integrated photoemission (PES) and x- 
ray absorption spectroscopy spectra (XAS) measurement. The experimental results are reproduced 
from Ref. [5] and the XAS results were shifted to obtain the best agreement with our theoretical 
results. 

Ml phase and the reduction of the band-width, we carried out LDA+DMFT calculations. 
Our theoretically calculated orbitally resolved spectral functions can be directly compared 
to the measured angle integrated photoemission (PES) and x-ray absorption spectroscopy 
spectra (XAS). In the last few years a large number of experimental studies of these kind 
were carried out probing also the many body character of the occupied and the unoccupied 
states in VO2 and serving as a stringent test of the theoretical approach. This validation 
of our LDA+DMFT results is crucial before proceeding to make reliable predictions for the 
strained materials for which these spectroscopic information is not yet available. 

The upfolded density of states of the LDA+DMFT, which includes besides the V-3d — t2g 
states the 0-2p and the V-3d — Cg states, are shown in Fig. I7j The dynamical correlation 
effects, not included in LDA, decrease the width of the quasiparticle t2g states for ~ 0.6 
times, in agreement with our calculated quasiparticle renormalization amplitude Z ^ 0.62. 
In the rutile structure two Hubbard bands appear: a very weak lower Hubbard band around 
-1.0 eV (was previously reported in theoretical [19j and experimental |S] studies) and a 
stronger upper hubbard band at 2.0 eV which overlaps with the Cg band (can be seen better 
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FIG. 8: (Color online) Upper panel: Brillouin zones and high symmetry points of the rutile (a) 
and the Mi (6) structures. Dashed arrows: the lattice vectors, thick arrows: the reciprocal lattice 
vectors, dashed lines: guide to the eye. Lower panel: The momentum resolved spectral function 
A{k, u) along the high symmetry paths (shown in the Brillouin zones by blue lines) of the R 
(left) and Mi (right) phases. For rutile the Brillouin zone corresponds to the doubled unit cell as 
indicated by the lattice vectors. The full lines show the LDA bands, and the color-coding shows 
the LDA+DMFT spectra, uj is in units of eV. 



in Fig. ^ where the Cg band is not shown). In the Mi phase, the most prominent effect of 
the DMFT theory is the appearance of a gap of ~ 0.7 eV at the chemical potential. For 
both the rutile and Mi phase the positions of the 0-2p and the V-3d — eg states show good 
agreement with the experiments. Our calculated Z factors for the Mi phase are ~ 0.9 for 
the e^ state and ~ 0.7 for the aig state confirming that the Mi phase is "less correlated" 
than the rutile phase, as emphasized in Ref. EH 

Fig.|8]shows the momentum resolved spectra in the R and Mi-phase. Notice the downshift 
of the two bonding bands in the Mi phase, primarily of x"^ — y"^ character (notice that there 
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are four V atoms per unit cell, and hence two bonding x^ — y'^ bands). The states above 
the Fermi level move slightly up and shrink due to many-body renormalization similarly to 
the t2g states in the R phase. The bands around 1 eV also acquire substantial lifetime in 
both phases. The oxygen bands below — 1.5eV and the eg bands above 2eV are almost 
unchanged compared to LDA. 

The results for the orbitally resolved 3(i — t2g spectral functions are shown in Fig. [9j In 
the rutile phase one can see that the three t2g states are approximately equally occupied 
predicting isotropic transport properties. For the aig state the bonding anti-bonding struc- 
ture can be also recognized similarly to the Mi phase with a splitting of ~ 0.4 eV which 
is by a factor of 2 smaller then in the LDA results. The calculated positions of the t2g 
states agree well with the XAS results. In the Mi phase, the weight redistribution is very 
different: only a single state, namely the bonding aig orbital is occupied. One can observe 
that the bonding aig state is shifted for ~ 0.8 eV lower, which is in good agreement with 
the previous theoretical results [12]. The spectral function at the upper edge of the gap 
has predominantly the e^ character, but there is also some weight of the anti-bonding aig 
character, which is in agreement with experiment showing that the spectral density has not 
purely e^ character above the gap [5J. The first two peaks at ~ 0.6 eV above the chemical 
potential are attributed mainly to the e^ states, and the third one is due to the aig state. 
This is consistent with the recent results of polarization dependent O K XAS experiments 
from Koethe at al. [3] where the orbital character of the states can be deduced by changing 
the polarization of the x-ray from parallel (O K XAS ||) to the c-axis to perpendicular polar- 
ization (O K XAS _L). The anti-bonding aig state lies at 1.3 eV above the Fermi level. This 
peak did not appear in previous theoretical studies of Ref fT9| but agrees well with current 
XAS results, which show that there is a prominent aig ~ 1.0 eV above the e^ peak [5|l39|H0]. 
It is interesting to note that the position of the anti-bonding aig state is roughly the same 
in both the DMFT results and the pure LDA results. The calculated separation between 
the bonding and anti-bonding peaks of the aig state is ~ 2.1 eV, which agrees reasonably 
with the experimentally found value (2.5 — 2.8 eV)|5]. Finally, the peak around 3eV of the 
XAS spectra may be assigned to the contribution from the Cg states, included in Fig. [7l 
but excluded in Fig. [o] This is consistent with the experimental finding of Ref. [5] where 
negligible change of the peak weight was observed across the MIT, but strong sensitivity to 
the polarization was noticed. 
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FIG. 9: (Color online) Comparison of the orbitally resolved t2g DOS of LDA+DMFT calculation 
corresponding to a V atom in rutile (left) and Mi (right) phases (r = 1.00) with angle integrated 
photoemission (PES) and x-ray absorption spectroscopy spectra (XAS) measurement. The experi- 
mental results are reproduced from Ref. [5j. The XAS result is shifted to obtain the best agreement 
with our theoretical results. 

To elucidate the effect of the strain along the c-axis the spectral functions of the t2g 



orbitals for different r ratios (r = 0.98, 1.00 and 1.02) are shown in Fig. 10 for both phases. In 
the rutile phase the el states are hardly affected by the strain. In contrast, the bonding anti- 
bonding splitting of the aig states shows a strong sensitivity to strain. The width of the upper 
Oig peak at ~ 0.5 eV decreases with increasing r, which is due to the weaker hybridization 
between the V atoms, as indicated by the decreasing hopping integrals. Surprisingly, the 
spectral weight at the chemical potential is practically unaffected by the changes of the 
lattice constant c. Similar calculation using cluster DMFT was carried out for the rutile 
phase, allowing for the formation of split bonding-anti-bonding pairs along the c-axis. We 
did not find an appreciable sign for the development of the bonding-antibonding splitting 
even for the r = 0.98 case, confirming that the single site DMFT is quite accurate in the 
rutile phase. 



Fig. 10 clearly shows that the width of the bonding aig peak is increased with decreasing 
r ratio, which can be attributed to the increase in inter dimer hoppings. For all ratios 
r, the gap in orbital space is indirect, i.e., the valence band is of aig character and the 
conduction band of el character. This is in agreement with the experimental findings of 
Ref. |5] demonstrating that the Peierls physics is playing a secondary role in the gap opening 
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in Ml phase. Due to the decreasing length of the c-axis, the e^ states are shifted to shghtly 
lower energy, which together with the broadening of the bonding aig peak results in the 
contraction of the gap, despite the increase in the bonding anti-bonding splitting of the aig 
peaks. The decrease of the gap size due to decreasing c-axis length is more apparent in the 



inset of Fig. 11 , which shows the total t2g DOS. This result is supported by the experimental 
result of Muraoka and Hiroi |9j demonstrating that the decrease in lattice parameter c leads 
to decrease in the metal-insulator transition temperature. This is a clear indication that 
smaller c-axis length leads to a weakened stability of the insulator in the Mi phase, and 
consequently a smaller gap in Mi phase. This behavior is not expected for a Peierls type 
gap, which increase as the lattice parameter decreases along the dimerized chains. 
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FIG. 10: (Color online) Orbitally resolved V— 3d — t2g DOS of a V atom in rutile (left) and Mi 
(right) phases in case of r = 0.98 (dotted, red line), 1.00 (full, black line) and 1.02 (dashed, green 
line). 



Fig. 



The evolution of the gap is reflected also in the gap of the optical conductivity |S] • In 
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the real part of the average optical conductivity aav = ^((^n + 2o"_l) is shown, where 
cr||(_L) is the optical conductivity in the case where the polarization of incident light is parallel 
(perpendicular) to the c-axis. The calculated optical gap and the intensity of the first peak 
corresponding to the t2g — t2g excitations compare well with the experimental results for 
polycrystalline VO2 films [H]. The shoulder around ~ 2.5 eV in the experimental optical 
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FIG. 11: (Color online) Real part of the averaged optical conductivity cjav = |(c|| +2o"_|_) of the Mi 
phase with different r ratios. Inset: total t2g density of states of the Mi phase. The experimental 
result is taken from Ref. 



conductivity, which is primarily due to the inter-band transitions, is shifted slightly upwards 
(~ 0.5 eV) in the theoretical result. This is an indication that the applied downfolding 
method describes well the low-energy properties, but not so well the higher energy interband 
excitations. 

The fact that the gap in the density of states arises between orbitals of different symmetry, 
indicates that the anisotropy of the transport properties in this phase will be very sensitive 
to disorder and grain boundaries which can drastically alter the orientation of these orbitals 
changing the matrix elements for hopping and conductivity. 



Fig. 12 shows the calculated and experimentally measured optical conductivity for dif- 
ferently polarized light. It can be seen that the trends of the dependence of the optical 
conductivity on the polarization are in a good agreement with experimental results, al- 
though the values are slightly different. When the polarization is perpendicular to the c-axis 
the optical response is practically unaffected by the strain. In case of parallel polarization 
the optical conductivity is strongly modified by changing the lattice parameter c, especially 
in the frequencies between 1.5 — 2.5 eV. This region can be attributed to the d^^-y^ — d^^-y^ 
excitations as can be concluded from the positions of the bonding and antibonding peaks in 



Fig. [TOj This results strongly indicate that the anisotropy of the transport properties is due 
to directed V-V bonds along the c-axis. 
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FIG. 12: (Color online) Real part of the optical conductivity corresponding to light polarized 
parallel ((Tii) and perpendicular (o"_l) to the c-axis in the Mi phase with different r ratios. The 
experimental result is taken from Ref. [14j (Exp. a) and \42\ I43j (Exp. b). 



In Fig. 13 , the orientation average optical conductivity of the R phase and the imaginary 
part of the self-energy at the Matsubara frequencies close to zero are shown. While the 
area below the calculated and experimentally measured [UJ optical conductivity (plasma 
frequency) agree fairly well, the width of the two Drude peaks are different. In order to im- 
prove the agreement with experimental results [H] an imaginary part of 0.55 eV (scattering 
rate) was added to the self-energy for the low frequency part of the optical conductivity, 
to simulate the experimentally measured broadening of the Drude peak. Inspecting the 



inset of Fig. [13} one can see that even by employing larger U values, the scattering rate 
(ImE(u; -> 0)) is not large enough to reproduce the experimental results, and the calculated 
optical conductivity will not show bad metal behaviour at this temperature. From this re- 
sult one can draw the conclusion that the experimentally measured large scattering rate is a 
consequence of an inhomogeneity of the system as indicated in recent experiments 
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FIG. 13: (Color online) Real part of the averaged optical conductivity aav = ^{(^W + 2(7_l) of the 
R phase. The experimental result is taken from Ref. [14j. Inset: Imaginary part of the self-energy 
corresponding to the t2g orbitals in the bonding anti-bonding basis at the imaginary Matsubara 
frequencies close to the real axis. 
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FIG. 14: (Color online) Sketch of the phase diagram of the rutile phase as a function of the 
temperature and the applied Coulomb repulsion U based on the calculated points (dots). In the 
green region (M) the metallic solution (shown in the left inset in terms of the Matsubara Green's 
function) is stable only while in the blue region only the insulator solution (right inset) can be 
found. In the red coexistence region (C) both solutions are stable. 
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IV. PHASE DIAGRAM AND LIMITS OF DOWNFOLDING 



In Fig. [T4| a sketch of the phase diagram of the rutile structure based on calculated 
points within the cluster DMFT is displayed as a function of the temperature and the 
Coulomb repulsion parameter U. Three different phases are distinguished: the metal (M), 
the insulator (I) and the coexistent (C) regions. To identify them, the imaginary part of the 
local Matsubara Green's function, linQ{un), of the cluster was investigated. We considered 
that the solution is insulating when lTn.Q{un) converges to at low imaginary frequencies for 



all orbitals, as shown in an example in the right inset of Fig. 14 The metallic solutions are 



the ones where Im^(u;„) tends to a finite value for at least one of the orbitals (see the left 



inset in Fig. 14). In the metallic region only the metallic solution is stable. In the insulator 
region only the insulating solution is stable, and in the coexistence region, both are stable. 
In order to decide whether a solution is stable or not, the DMFT calculations were started 
from an ansatz of a specified type (metal or insulator) and if it remained of the same type 
after the self-consistent solution is reached, one can regard it as a stable mean-field solution. 



In Fig. 14 it can be observed that below U ~ 2.9 eV only the metallic solutions are stable 
while above ~ 3.8 eV only the insulator is stable. The nature of the insulating state in 
rutile phase is very different from the insulating state in Mi phase. The lower Hubbard 
band in the rutile insulator is an almost equal mixture of all three t2g orbitals. However, 
the interaction strength needed to open the true Mott gap without the help of the Peierls 
mechanism, is considerably larger. In the coexistence region, one can expect a crossover 
between the two phases governed by the free energy of the system ^5]. There is a strong 
experimental indication that the rutile phase resides in the vicinity (but on the metallic side) 
of this crossover. Pouget and Launois showed that the metallic feature of the rutile phase is 
very sensitive by substitutional alloying of VO2 with Nb (Vi_2;Nba;02) which increases the 
c/a ratio and results in the appearance of a gap keeping the rutile structure at x = 0.2 |^ . 
Recently, Holman et al. |17] reported insulator to metal transition in Vi_^Moa;02 system 
at X ~ 0.2. All those experiments suggest that VO2 is in the crossover region near the 
coexistence of two solutions in the cluster DMFT phase diagram. 

On the basis of our calculation we conjecture that the rutile phase might be able to 
support either metallic or insulating solution (with a very small gap), and hence either of 
the two phases can be stabilized depending on small external stimuli. On the other hand. 
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the Ml phase at the same interaction strength supports only insulating solution, and it is 
unlikely that small external perturbation can turn it to metallic state. 

One possibility to improve the agreement between the experimental and theoretical results 
is that one should use different parameters in the downfolded model. The second possibility 
is that the material is strongly inhomogeneous which was not taken into account so far 
in any theoretical calculations. Finally, it is most likely that calculation without invoking 
the downfolding approximation will result in a more accurate description of this material. 
This would place the VO2 closer to the Mott charge transfer insulator boundary in the 
Zaanen-Sawatzky-Allen phase diagram. Work in this direction is in progress. 

V. CONCLUSIONS 

Our exploratory theoretical research set up the machinery for describing the subtle inter- 
play of Coulomb correlations, orbital degeneracy and strain in determining the mechanism 
of the MIT in VO2. Our theory, coupled with existing strain experiments, clearly shows 
that the Peierls distortion is only one element affecting the MIT and the switching mecha- 
nism of this material. The LDA+DMFT calculations in the unstrained material are in good 
agreement with experiments. We performed the first LDA+DMFT studies of the electronic 
structure of VO2 under strain. Besides the Peierls increase in aig bonding-antibonding split- 
ting, the lowering in energy of the e^ orbital, and the rapid change in bandwidth of the aig 
orbital due to the varying overlaps, play an equally important role in controlling the posi- 
tion of the MIT. These theoretical insights can be used for understanding and improving 
material properties by means of chemical substitutions. For a more accurate description it 
is mandatory to take into account the oxygen degrees of freedom, and calculations of the 
total energy in the regions suggested by this exploratory work. 
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